Lattice QCD Constraints on Hybrid and Quark Stars 
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A QCD-motivated dynamical-quasiparticle model with parameters adjusted to reproduce the 
lattice-QCD equation of state is extrapolated from region of high temperatures and moderate bary- 
onic densities to the domain of high baryonic densities and zero temperature. The resulting equation 
of state matched with realistic hadronic equations of state predicts a phase transition into the quark 
phase at higher densities than those reachable in neutron star interiors. This excludes the possibility 
of the existence of hybrid (hadron-quark) stars. Pure quark stars are possible and have low masses, 
small radii and very high central densities. Similar results are obtained for a simple bag model 
with massive quarks, fitted to reproduce the same lattice results. Self-bound quark matter is also 
excluded within these models. Uncertainties in the present extrapolation are discussed. Comparison 
with standard bag models is made. 
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I. INTRODUCTION 

Nowadays it is commonly accepted that the quark- 
gluon phase of the matter can be formed in the course 
of heavy-ion collisions at ultrarelativistic energies. There 
are indirect signatures that this state has been produced 
at Relativistic Heavy-Ion Collider (RHIC) 0. On the 
other hand, it was suggested long ago that quark mat- 
ter may exist either in the interior of some stars or as a 
new family of pure quark stars. Then it was realized that 
some rather massive neutron stars may have a quark core 
(so called hybrid stars) , cf . Refs d 0] and references 
therein. Recently the idea of quark pairing with high 
values of the gap (up to A ~ 100 MeV) attracted much 
attention, cf. Rcf. 5]. Observable signals of quark mat- 
ter, presumably existing in interiors of hybrid and pure 
quark stars, are under extensive discussion. The most 
promising of these signals seem to be an abnormal mass- 
radius relation [(|, r-mode spin-down and a specific 
steller cooling properties 

Although there is no a priori obstacle to existence of 
the quark matter in above the mentioned forms, the ac- 
tual theoretical predictions are very uncertain. They 
hinge on the lack of understanding of the confinement 
mechanism in QCD. Therefore, parameters of theoret- 
ical quark models cannot be anchored to experimental 
observables. Most of the treatments are based on ei- 
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ther the standard or slightly modified MIT bag model, or 
the Nambu-Jona-Lasino (NJL) model, cf. Ref. and 
references therein. In terms of a bag model, the main 
uncertainty is associated with the value of the vacuum 
pressure, i.e. the bag constant B c g. The steller proper- 
ties are sensitive to this. A large value of B g excludes 
the existence of hybrid stars. With a decrease of J3 e ff, 
first the most massive neutron stars acquire a quark core, 
becoming hybrid stars. For a small -B c ff, strange quark 
matter becomes absolutely stable, i.e. more stable than 
56 Fe [Ulnl]. Then pure quark stars (strangelets) of an 
arbitrary size may exist, being hold together by strong 
rather than gravitational forces. Physical implications of 
the existence of a self-bound strange quark matter are 
reviewed in Ref. [l2^ . 

A straightforward calculation of the properties of 
strongly interacting matter and, in particular its equa- 
tion of state (EoS), is possible on the lattice 13]. This 
technique is progressing rapidly. So far, the main body 
of lattice results concern the case of zero baryon chemi- 
cal potential (/i = 0) and finite temperatures. This case 
is just opposite to that of cold baryon matter present in 
the quark stars. Recently new lattice predictions on the 
EoS at finite, but small, baryon chemical potentials be- 
came available 

El 13 El Nevertheless, they are still 
far from the astrophysically relevant domain. Therefore, 
a theoretical modeling is needed for their extrapolation. 

The interpretation of these lattice results within phe- 
nomenological quasiparticle models, i.e. in terms of effec- 
tively massive quarks and gluons with a simple interac- 
tion, turned out to be very successful both at zero chemi- 
cal potential [IIEIEIIMMJ2IJIJMIMIMI2II2I 
and at small finite ones |25l l2q. 1271 128| . With few phe- 
nomenological parameters it was possible to reasonably 
reproduce all lattice thermodynamic quantities. Effec- 
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tively the results can be interpreted in terms of a modi- 
fied MIT bag model for massive quarks with density and 
temperature dependent bag constant and masses. These 
models can be used to extrapolate the lattice EoS to the 
domain of cold baryon matter. This has been done by 
Peshier, Kampfer and Soff [211 for the thermodynamic 
quasiparticle model of Ref. [2(| ■ They considered the pos- 
sibility of the existence of quark stars and found a mass- 
radius relation for them. They concluded that quark 
stars with a mass larger than 1M© may only exist for val- 
ues of the effective bag constant B 1 / 4 < 180 - 200 MeV, 
whereas their model produced a higher value. The pos- 
sible existence of hybrid stars was not considered in that 
work. It is important to mention that the quasiparticle 
model of Ref. [2l| was fitted only to the fj, = lattice 
data available at that time and then it was extrapolated 
to finite values of the baryon chemical potential. The 
later fit of the same model to recent lattice data at finite 
baryon chemical potentials [14j resulted in somewhat dif- 
ferent parameters of this quasiparticle model (25|. The 
same quasiparticle model of Ref. [2(J was used in Ref. [24[ 
in order to extrapolate the HTL calculations to the cold- 
baryon-matter domain. Conclusions very similar to those 
in Ref. |2l| were derived concerning the properties of 
quark stars. 

In this paper we would like to make "lattice QCD mo- 
tivated" predictions on the possible existence of hybrid 
and quark stars, as well as their properties. To do this, 
we extrapolate the existing lattice results to the region 
of cold electroneutral baryon matter relying on the phe- 
nomenological QCD motivated dynamical-quasiparticle 
(DQ) model recently proposed in Ref. |2j|. In equi- 
librium, at high temperatures the DQ model complies 
with the quas iparticle picture of the hard thermal loop 
approach 29], whereas at lower temperatures it simu- 
lates the confinement of the QCD. The latter is mani- 
fested by the fact that the solution to the model equa- 
tions simply does not exist below certain combination of 
the temperature and the chemical potential. Two sets 
of parameters of this model (below denoted as "1-loop" 
and "2-loop") were fitted to reproduce lattice data at 
finite baryon chemical potentials Extrapolation of 

the lattice data to the case of the cold baryon matter is 
essentially simpler in this model than in that of Ref. [2lJ , 
where an implicit thermodynamic-consistency equation 
should be numerically solved. In fact, this advantage 
follows from an explicitly thermodynamically consistent 
formulation in terms of dynamic rather than thermody- 
namic variables. 

In Sect.[n]we start with brief recapitulation of the DQ 
model. To have a reference point for the DQ model, we 
also consider three versions of the MIT bag model. Two 
"light-bag" models ("light-bag-155" and "light-bag-200" ) 
with conventional current quark masses (m„o = 5 MeV, 
mdo = 7 MeV and m s0 = 150 MeV), zero gluon mass 
(m g = 0) and different bag constants (B 1 / 4 = 155 MeV 
and 200 MeV) covering a broad range of EoS 's usually ap- 
plied to the treatment of hybrid and quark stars |3(J, |31| . 



We also consider a "heavy-bag" model, the quark/gluon 
masses and bag constant of which are fitted to reason- 
ably reproduce the lattice results of Ref. 0]: TO « = 
m d = 330 MeV, m s = 450 MeV, m g = 600 MeV and 
5 1 / 4 = 183 MeV. In Sect.lTTTl we confront the predictions 
of all these models to the (2+l)-flavor lattice data 0] 
at finite chemical potentials. In Sect. HVl we derive the 
predictions of all the above models regarding the possi- 
ble existence of hybrid and quark stars and discuss their 
properties. 



II. MODELS 
A. Dynamical-Quasiparticle Model 

The effective Lagrangian of the DQ model [2^| treats 
transverse gluons </> a and quarks ?/>/ c of Nf = 3 flavors 
and N c — 3 colors. For this model it is essential that 
Nf = N c . These constituents interact via mean fields 

C = {Cu, Cd, Cs}' 

L 0=1 
N c N f 

+ E E tf&i^ - m f(CMfc - b(x), (i) 

c=l f=l 

where N g = 2(N 2 — 1) is the number of transverse gluons, 
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Here m/o is the current /-quark mass, whereas m g and 
rrif are effective masses of gluons (g) and quarks of flavor 
/, respectively. The masses and the coupling constant 
g 2 (x) depend on mean fields 



B(x) = b c - j^r [xV (x) - xhg 2 (xc)] 

2 f x 

+ ~rr / d X ix\g 2 {xi) , 



(•5) 



is the potential of the mean-field self-interaction having 
a meaning of an effective bag parameter. Here Be and 
Xc are parameters of the model. 

The equations of motion for the mean fields read 

dC dB 1 dm 2 Q 2 l ^ dml 



3 



where 



N„ 



(ipfcfpfc) 



a=l 

* — ' rat 

c=l 3 

N r r°° k 2 dk 



2,t j ./ (fc 2 + m2)i/2 /9(fc) ' (7) 



n2 J o (fc9 + ro , )i75 [/,./(*) +/,./(*)]. W 

are scalar densities of gluons and quarks (divided by the 
mass), respectively, and f g (k), f q j(k) and f q j(k) are 
occupation numbers (distribution functions) of gluons, 
quarks and antiquarks. In the case of equilibrium that 
we consider here, these are 



fM = 

/*,/(*) = 
f q ,f(k) = 



1 



eM(k 2 + m 2 g ) 1/2 /T]-l' 
1 

exp{[(fc 2 + to 2 ) 1 / 2 - fi f ]/T} + 1 ' 
1 

exp{[(fc 2 +m2)i/2 +M/ ]/T} + l- 



(9) 
(10) 

(11) 



where T is the temperature, and fif is the /-flavor quark 
chemical potential. In the general case, all /if may be 
different. 

The solution of above equations of motion © with 
respect to the mean fields is 



= rf 



AN, 



(12) 



Here we heavilyrely on the condition N c = Nf. It is 
easy to verify |2g that with this solution the effective 
masses of gluons and quarks, cf. Eqs 01 and |0J, re- 
produce the hard-thermal-loop results |29| in the high- 
temperature limit, provided the coupling constant g 2 (x) 
is appropriately defined. In fact, this property was the 
main requirement in construction of this DQ mod el l28| . 
The appropriate choice of the coupling constant [28j is 
the following 



9 2 (X) = 



16tt 2 



PoH(.x 2 + xt)/x 2 



J fix) , 



(13) 



where /3 = |(11JV" C - 2Nf). Nf is the effective num- 
ber of quark flavors at the energy scale x, which may 
be Nf < Nf ;65]. The parameter Xc ^ s introduced in 



Eq. (J5J and Xo is another phenomenological parameter 
of the model. An auxiliary function /(x), satisfying the 
condition f(x ~^ °°) ~* 1; helps us to choose between 
the 1-loop and 2-loop perturbative limits of the coupling 
constant |32ll33|. 

The "1-loop" version of the model 

/l-loop = 1 (14) 
overestimates by approximately 10% the lattice data (see 
Ref. [H and discussion in Sect. ITTTj) . Therefore, the "2- 
loop" version was also considered. The latter was chosen 
in the form 



/2-ioo P (x) = 1 + arctan 



in 2 P 



! (x)ln 



9 2 (X) 



(15) 



with 
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13iV c 7V| ff 



■ 3Nf/N c 



A = 0.001 



167T 2 



With this /2-ioop, the coupling constant meets the 2-loop 
perturbative limit at x —> 00 ■ The additional g 2 ln(0.001) 
term is sub- leading as compared to g 2 (x) In g 2 (x) and 
thus does not prevent agreement with the 2-loop approx- 
imation for the coupling constant. In fact, the function 
/2-ioop is an "exotic" representation of a constant, since 
in the range, covered by lattice simulations, /2-ioo P (x) — 
2.6 with good accuracy. Precisely this enhancement of 
the coupling constant is required to fit the actual overall 
normalization of the lattice results. 

As we have mentioned the present model simulates the 
confinement of quarks and gluons. At some value of x 
the argument of ln[(x 2 + Xo)/Xc\ m ^1- becomes 
very close to 1, and hence g 2 — > 00. Due to that there 
are no solutions to the above equations below certain 
values of temperature and chemical potentials. This can 
be interpreted as a kind of confinement. 

To summarize, the procedure of solving the model 
equations is as follows. First, we define all the free pa- 
rameters of the model (xc, Xo, Be), including the aux- 
iliary function /(x)- Given the temperature T and the 
set of chemical potentials /i / , the set of equations (J2J) © 
and Q~ Q«3) should be solved. As a result of this solu- 
tion, we obtain effective quark and gluon masses and the 
value of x variable, which is required for calculation of 
B(x), cf. Eq. JSJ. When all the quantities are defined, 
we can calculate the energy density e(T, {/i/}), pressure 
P(T, {/df}), and baryon density n(T, {^/}) as follows 



E (T,{A*/}) = (l-^^J^k 2 dk(k 2 + m 2 g )'/ 2 f g (k) 

Nf at roc 

11 Jo 



(16) 
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Version 


T c a , MeV 


m u , MeV 


m d , MeV 


m s , MeV 


m g , MeV 


B 1/4 , MeV 


c 


light-bag- 155 


167 


5 


7 


150 





155 


0-^0.3 


light-bag-200 


167 


5 


7 


150 





200 





heavy-bag 


167 


330 


330 


450 


600 


183 


0^0.3 



"The critical temperature Tc is chosen on the condition of the 
best reproduction of lattice predictions 1 1 -41 by the "heavy-bag" 
model. For the "light-bag" models it was simply kept the same. 



TABLE I: Parameters of bag models 
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JV f 



(17) 
(18) 



The gluon and antiquark terms vanish at zero temper- 
ature. Here we have introduced factors (l — |c) and 
(1 — c) to make provision for the bag model (see the 
next Subsect.). For the DQ model, c is identically zero, 

CDQ-modcl = 0. 

In Eqs 1)16(1 and 1(17(1 we have also made provision for 
a pairing contribution which may be important at 

low temperatures. The respective term can be directly 
added to the Lagrangian. It preserves the above mean- 
field solution, provided it is independent of these mean 
fields. For temperature higher than a certain temper- 
ature T c csc , this contribution vanishes, SB a — 0. For 
T < T c csc the quarks can be paired, forming a color su- 
perconductor. This gives rise to SBa that depends on 
the color superconducting gap. The form of this term, 
SB a, depends on the pairing channel. There exists a 
variety of possible phases. We will present the results 
for two phases usually discussed. One expects that the 
phase that can be realized at moderate densities is the 
two-color-superconducting (2SC) phase. At larger densi- 
ties the most symmetric color-flavor-locked (CFL) phase 
may be formed (for A > 3m 2 /(2/i)). In these cases we 
have 

SB A ~ - l/A J* at T < T c csc , (19) 

where A(T) is the gap, v — 1/9 for 2SC phase and 
v = 1/3 for CFL phase, and [i is the baryon chemical 
potential. 

Note that thermodynamic consistency is automatically 
fulfilled in our scheme, since we proceed from a proper 
Lagrangian formulation. The presented model is of pure 
mean-field nature. In particular, it disregards fluctua- 
tions that may be important in the vicinity of both the 
deconfinement and the superconducting phase transitions 



35]. However, these fluctuations are usually disregarded 
in other models applied to both high and low tempera- 
tures as well. These effects require further study. 



B. Bag models 

The last formulas of the previous subsection, i.e. Eqs 
(II 01) 1(18(1. are directly applicable to the bag model with 
the only simplification that all masses m u ,m,d,m s ,m g 
and the bag constant B are really constants, i.e. are in- 
dependent of any mean fields, temperature and chemical 
potentials. Often perturbative corrections oc q 2 to the ki- 
netic terms of both quarks [3(| and gluons j3]} given by 
the factor c (which is 2a 2 /tt in case of massless quarks 
and gluons) are also taken into account, see Eqs 1(16(1 - 
((18(1 . Since perturbation theory is not applicable to the 
whole region of temperatures and chemical potentials of 
our interest, we will not use an explicit expression for this 
correction but rather vary it within a reasonable range, 
c = -=- 0.3, following Ref. j^. 

Below we consider three versions of the bag model, the 
parameters are listed in Table [I] The two first versions, 
labeled "light-bag-155" (similar to the model used in Ref. 
[3(j) and "light-bag-200" (similar to the model used in 
Ref. represent a range of "conventional" bag mod- 

els used in astrophysical applications. As we will show 
below these models cannot reproduce lattice results. The 
"heavy-bag", as well as DQ based models, are fitted to 
reasonably reproduce lattice data of Ref. [lij. see next 
section. This fit required heavy quarks and gluons to be 
introduced. 
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"1 loop" 

"2 loop" 
heavy-bag 
ght-bag-200" 
light-bag-155" 



0. 

CO 



T/T n 




FIG. 1: (Color online) Pressure (left panel) and the interaction measure, e — 3P, (right panel) scaled by T 4 as functions of 
T/Tc at zero baryon chemical potential /i = 0. The solid line corresponds to "1-loop" version, whereas the dashed line, to 
the "2-loop" one of the DQ model. The vertically hatched band is the "heavy-bag" fit with c = and c = 0.3 for the top 
and bottom boundary lines, respectively. The grey band (for P) and the grey line (for e — 3P) represent prediction of the 
"light-bag-155" model. The dash-dotted line corresponds to the "light-bag-200" model. Lattice data, triangles and dots, are 
from Ref. 0]. 



III. COMPARISON WITH LATTICE RESULTS 

In this section we demonstrate the quality of the fit 
of the above described models to the recent (2+lVflavor 
lattice predictions at various chemical potentials 

We again start with a brief recapitulation of Ref. [28j , 
where this comparison was done for the DQ model. To 
be consistent with the lattice data, current quark masses 
m-no = m-do = 65 MeV and m s o = 135 MeV were adopted, 
as in these lattice calculations. It was found that the 
actual results of our quasiparticle model are quite insen- 
sitive to variation of rrifo from above lattice values to 
the "physical" ones m U Q = to^o = 7 MeV and m s o — 
150 MeV. This insensitivity is to the extent that predic- 
tions of the DQ model with lattice and "physical" quark 
masses are hardly distinguishable in the scale of figures 
presented below. In fact, this is not surprising, since 
the 80-^90% of the effective quark masses (0J result from 
interaction with the mean field B(\)- The model in- 
volves several phenomenological parameters: the param- 
eter Be, cf. (J5J, the "QCD scale" xc, and an auxiliary 
function f(x), cf. (|T3|) . Another parameter Xo should be 
taken small Xo ^ Xc- ^ shifts the lower limit of inte- 
gration in the expression for B(x), cf. Eq. JSJ, from the 
singular point of the coupling constant, cf. Eq. l(T5jl. and 
hence regularizes the calculation of B(x)- Therefore, it 
is closely related to the parameter Be, which is an inte- 
gration constant in the same expression. A change of Xo 
implies the corresponding change of Be- In all further 
calculations we take Xo = ^■^Xc> an< ^ the va hies of Be 
stated below correspond only to this choice. 

An implicit parameter of the model is the critical tem- 
perature Tc, i.e. the temperature at which the decon- 
finement phase transition occurs at fi = 0. It would not 
be reasonable to identify this temperature with that of 



the end point of the solution discussed above. The rea- 
son is that the phase transition at /i = in the case of 
(2+1) flavors is of the cross-over type, as found in lat- 
tice calculations. This implies that a strong interplay 
between quark-gluon and hadronic degrees of freedom 
occurs near Tc, which actually determines the Tc value 
itself. Since hadronic degrees of freedom are completely 
disregarded in the model, it cannot properly determine 
the Tc value. Therefore, the value of Tc was varied from 
the determined end-point temperature to slightly below 
in order to achieve the best fit of the lattice results. The 
fitted value Tc — 195 MeV is slightly above its lattice 
counterpart 175 MeV. 

As demonstrated in Ref. [28|, the "1-loop" version of 
the DQ model, cf. Eq. i|14|) . does not reproduce the over- 
all normalization of the lattice pressure. For the best 
fit of the lattice data the overall normalization factor 
(-Vnorm) was chosen equal 0.9 for the "1-loop" calcula- 
tions. The value -/V norm is indeed somewhat uncertain. 
Lattice calculations were done on lattices with N t = 4 
temporal extension [l4|. To transform the "raw" lattice 
data into physical values, i.e. to extrapolate to the con- 
tinuum case of Nt —>■ oo, the raw data are multiplied 
by "the dominant T — > oo correction factors between the 
N t = 4 an d continuum case", c p = 0.518 and c M = 0.446, 
[l4| . These factors are determined as ratios of the Stefan- 
Boltzmann pressure at fi = (c p ) and the ^-dependent 
part of the Stefan-Boltzmann pressure (c„ ) to the cor- 
responding values on the N t = A lattice [14(- In view of 
this uncertainty, it is reasonable to apply an additional 
overall normalization factor ./ Vnnrr n to the same quantities 
calculated within DQ model |66|. In order to keep the 
number of fitting parameters as small as possible, a sin- 
gle normalization factor A^orm was applied instead of two 
different ones, c p and c M , as in the lattice calculations. 
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Version 


Nf 


T c , MeV xc, MeV Bc/xc" 


/-factor 6 TVnorm 


"1-loop" 


3 


195 141.3 -97.5 


1 0.9 


"2-loop" 


3 


195 119.6 -267.5 


2.6 1 



"These Be values correspond to the choice Xo = O.Olx^- 
'This is the effective value of the auxiliary function f(x) in the 
temperature range under investigation, i.e. from Tc to 3Tq- 



TABLE II: Best fits of DQ-model parameters to lattice 
data [3 



As for the "2-loop" version, the auxiliary function 1|15|) 
was chosen to reproduce the lattice data normalization. 
In fact, the function /a-ioop is an "exotic" representation 
of a constant, since in the temperature range under con- 
sideration, from Tc to 32c, /2-ioo P (x) — 2.6 with good 
accuracy. In this respect, any function /, providing us 
with the additional factor 2.6 in the temperature range 
from Tc to 3Tc and respecting the proper perturbative 
limit of the coupling constant at very high temperatures, 
is just as suitable for this fit. The two sets of parameters 
of the DQ model are summarized in Table ITTI 

The comparison with (2+l)-flavor lattice results 0] 
of the above two versions of the DQ model, as well as 
the various bag models of Tabled is presented in Figs.Q] 
and |21 The baryon chemical potential fi is defined as 
follows: fj, u = Hd = A*/3 and ji a — 0. Both versions of the 
DQ model perfectly fit the lattice data above Tc- This 
fact demonstrates that the way to fit these data even 
within the DQ model is not unique. Moreover, as it was 
demonstrated in Ref. [2^, even 2-flavor lattice pressure 
(l5| is well reproduced above Tc by the DQ model with 
practically the same parameters. 

For two of the bag models we have taken into account 
the uncertainty associated with the perturbative correc- 
tion c. Following Ref. [3(J we varied c from to 0.3, see 
Table [I] Therefore, the results for bag models ( "light - 
bag-155" and "heavy-bag") are displayed, as a rule, by 
bands with c = and c = 0.3 for the top and bottom 
boundary lines, respectively. In the case of the "interac- 
tion measure" e — 3P and light quarks, the kinetic terms 
cancel with good accuracy and the band shrinks to a 
line. For the "light-bag-200" model, we do not display 
uncertainties associated with c in order not to overload 
the figure and since the general picture seems clear with- 
out this additional variation: the "light-bag" models do 
not follow lattice predictions. We checked that these fits 
to the pressure and "interaction measure" e — 3P could 
be improved at the price of a significant increase of the 
bag constant B. However, reproduction of these lattice 
quantities still remains rather poor even with this im- 
provement. Moreover, the remaining quantities, AP and 
riB, prove to be unchanged, since they are independent 
of B. Within the above mentioned c-uncertainty, the 
"heavy-bag" model reasonably well covers the lattice data, 
however, at expense of introducing large quark and gluon 
constituent masses. 



To reveal similarities between the DQ and bag models, 
in Fig. we compare quark/gluon masses and bag pa- 
rameters from these two kinds of models. Note that these 
are the only quantities that matter in the calculations of 
thermodynamic functions such as those of Eqs (|16|) - (|18f) . 
The u and d quark masses of the DQ model are indis- 
tinguishable within the scale of this figure. Quark and 
gluon masses, required for the reproduction of the lattice 
results, are high (> 300 MeV) in both the DQ model and 
the "heavy-bag" model. Thus, the failure of the "light- 
bag" models to reproduce lattice data can be associated 
with the light quark and gluon masses used in the calcu- 
lation. Note that in the DQ model, B becomes negative 
at high temperatures. This means that the interaction 
changes from repulsion to attraction. This behavior of 
the bag parameter is common for models reproducing 
lattice data, cf. Refs El23]|. 

As argued in Ref. [38|, the lattice data below Tc are 
well reproduced by the resonance hadronic gas model. 
However to comply with the lattice results, hadron 
masses in this resonance hadronic gas should be en- 
hanced. In particular, the pion mass was taken ~ 
770 MeV, as resulted from lattice calculations. This fact 
indicates that the EoS in the hadronic sector is not reli- 
ably predicted by the lattice calculations. The analysis 
of Ref. [3!| shows that strong interaction effects might 
be crucially important in the description of the hadron 
system at high temperatures and small baryon chemical 
potentials. Also we cannot use the resonance hadronic 
gas model at the low temperatures of present interest 
due to the lack of interactions incorporated in this model. 
This is why we avoid fitting our realistic hadronic EoS's, 
discussed in the next section, to the lattice data. 



IV. EOS AT T = 0: HYBRID AND QUARK 
STARS 

Possible existence of compact stars consisting either in 
part (hybrid stars) or completely (quark stars) of a decon- 
fined quark matter is one of the vividly discussed asp ects 
of modern astrophysics, see recent reviews [iol l4ll 142^. 
Typically, for the description of a deconfined phase a cer- 
tain version of the MIT bag model is used, see Ref. [43] for 
a survey of parameterizations. In this section we explore 
a possibility for hybrid and quark stars, proceeding from 
the DQ and heavy-bag models for the EoS of deconfined 
matter. 



A. Cold star matter 

For matter in /3-equilibrium at zero temperature sev- 
eral conditions should be met. Chemical potentials of 
quarks and leptons satisfy equilibrium conditions with 
respect to reactions u + e <-» s, d <-» s, /1 <-» e: 

Mu + Me = Ms , Md = Ms . Mm = Me ■ (20) 
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FIG. 2: (Color online) AP = P(T, fi) - P(T, fj, = 0) (left panel) and the baryon density n (right panel) scaled by T 4 and T 3 
as functions of T/T c at nonzero baryon chemical potentials fj, = 330 and 530 MeV (from bottom to top). Note different scales 
(left and right, respectively) for /i =330 and 530 MeV in the right panel. Notation is the same as in Fig. 
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The baryon (neutron) chemical potential is related to the 
quark chemical potentials as 



(21) 



In Ref. it was suggested that in multi-component 
systems the first-order phase transition proceeds through 
a mixed phase. In our case this would imply a coexistence 
of quark and hadronic matter in a certain density range. 
In the mixed phase both quark and hadron subsystems 
are charged and electro-neutrality is realized only glob- 
ally (see also Refs [HEH)- Coulomb and surface effects 
were disregarded in those works. The latter effects, how- 
ever, play a crucial role. As argued in Ref. 0], the re- 
gion of the hadron-quark mixed phase is very narrow, if 
not absent at all, due to rather high surface-tension and 
charge-screening effects on the quark-hadron boundary. 
Therefore, we can use Maxwell construction in order to 



describe the quark-hadron phase transition 67] . This im- 
plies the local electro-neutrality 



1 



1 



3 n - - r d - r s 



n e ~n l _ l = 0, 



(22) 



where n u , n^, n s are quark densities, n e is the electron 
density and is the muon density. The lepton contri- 
bution to the pressure is given by 



riifj,i 



p 2 dpy / p 2 + m\ 



(23) 



where m e and m M are lepton masses, jUj = y/ppA + m 2 , e 
and n densities are 



(li 2 - TO 2 ) 3 / 2 

[jh 90m -rm), 



(24) 
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and 6(x] is the ordinary step function. 

Now let us specify the hadronic EoS. The Urbana- 
Argonne EoS ||g ,418 + Sv + UIX* complies with the 
rich experimental information available for n < uq (uq 
is the nuclear saturation density). For that reason this 
EoS is considered to be the most realistic EoS at present. 
However, it uses a non-relativistic potential and hence 
violates causality at n > 4n . Ref. [3l| suggested a 
simple analytical parameterization (below referred as the 
HHJ EoS) of A18 + 5v + UIX* EoS valid for n < 4n , 
which respects causality at higher densities. On the other 
hand, as an extrapolation to higher densities, relativis- 
tic mean field (RMF) models are practical. However all 
RMF models have an unpleasant feature. They produce 
a large fraction of protons in neutron star matter that 
permits a very efficient cooling process, i.e. the direct 
Urea (DU) process n — » peD, for rather low densities 
< 3no- First, this disagrees with predictions of the men- 
tioned A18 + Sv + UIX* EoS, where the DU process is 
allowed only for n > 5tiq. Second, a low threshold den- 
sity $3«o seems to be at odds with a neutron star 
cooling scenario 0] predicted by the Urbana-Argonne 
EoS. 

Another problem is that rather massive neutron stars 
may exist. Recent measurements of neutron star masses 
in bina ry c ompact systems yielded Mjo75i+i8Q7 = 2.2 ± 
O.2M |5(j at the 95% confidence level. Even assuming 
the highest value of the uncertainty of this measurement 
one can conclude that MJ0751+1807 > 1.6M@ (within 2a). 
Some models with a so-called soft hadron EoS do not 
satisfy this requirement. Note that the Urbana-Argonne 
EoS yields the limiting mass of a neutron star ~ 2.2A/©. 
With the relativistic HHJ EoS the limiting neutron star 
mass decreases slightly below ~ 2 M©. Possible phase 
transitions to the cither pion or kaon condensation, hy- 
peronization, quark matter, etc. soften the EoS, that 
may significantly decrease the limiting mass of the star, 
cf. Ref. 0. 

Ref. fUj suggested a solution of both problems within 
an effective RMF model with field dependent effective 
hadron masses and coupling constants. We will use two 
models from this work. One model, "hadron-MW2" 
(MW(nu)(z = 2.6) in the notation of Ref. [13), fits the 
HHJ EoS for all densities yielding the threshold density 
for the DU process ~ 5.2uq. Another model, "hadron- 
MW1" (MW(nu)(z = 0.65) in the notation of Ref. H3) 
fits the HHJ EoS for n < 2no but allows to increase the 
limiting mass of a neutron star due to a stiffening of the 
EoS at a higher density. The threshold of the DU pro- 
cess is ~ 4no which does not contradict to the above 
mentioned description of the star cooling. Thus, the 
"hadron-MWl" model possesses the same advantages, as 
the Urbana-Argonne EoS and the "hadron-MW2" one, 
but in addition allows to increase the limiting mass of 
the star. The pressure of the "hadron-MWl" model as 
a function of the baryon chemical potential is lower as 
compared to that of the "hadron-MW2" model, which is 
in favor of the phase transition to the quark state at a 



smaller density. In order to keep the considerations sim- 
ple, we use models that disregard the possibility of p~ 
condensation discussed in Ref. |H^ . 

In Fig. 0] (left panel) we show effective masses of quarks 
at T — as a function of the baryon density in units 
of the normal nuclear matter density tiq ~ 0.16 fm -3 . 
Horizontal lines represent quark masses of the "heavy- 
bag" model, see Table [I] Self-consistent calculations of 
the DQ model produce a much smaller difference between 
strange and non-strange quark masses compared to 120 4- 
150 McV, which is usually accepted by bag models. In the 
DQ model, quark masses slightly decrease with increasing 
density up to (8-9-9) no- For higher densities the behavior 
reverses and the masses undergo a slight increase. Note 
that such a density behavior of the quark masses is quite 
different from that obtained in the NJL model [53|]. In 
the latter case for zero temperature the quark masses 
drop with increasing density, reaching zero at a moderate 
value of the baryon density n c , they also drop to zero at 
T = Tc- Therefore it seems doubtful that the NJL model 
could reproduce the lattice data at least near and above 
Tc\ similarly we have shown that the "light-bag" models 
fail. 

Fig. El (right panel) demonstrates the behavior of the 
effective bag parameter B 1 / 4 as a function of density at 
zero temperature. The B values prove to be very high in 
the DQ model, whereas the bag parameter of the "heavy- 
bag" model is in the range of ordinarily used values. Im- 
plications of a temperature-density dependent bag pa- 
rameter, in particular for quark stars, have been earlier 
considered in Ref. 0, 13] . 



B. Do hybrid stars exist? 

In Fig. we present the pressure as a function of the 
baryon chemical potential (left panel) and as a function 
of the baryon density (right panel) for the models con- 
sidered above. Realistic hadronic RMF EoS's are repre- 
sented by the "hadron-MW2" and "hadron-MWl" mod- 
els discussed above. The latter EoS is stiffer, as is seen 
from the right panel. As a reference model, the stiffest 
RMF EoS of the original Walecka model ( "hadron- W" ) is 
also demonstrated. The crossing points of the hadronic 
and quark curves on the (P, /1) plane (left panel) show 
critical values (/i c , P c ) for occurrence of the first-order 
phase transition from the hadronic to the quark phase. 
The latter transition arises, if for fi > fj, c the pressure 
of the quark phase is higher than that for the hadron 
phase. The corresponding critical densities and n q c 
at the Maxwell construction, n£ < n% are listed in Ta- 
ble EH 

As seen, all latticc-QCD-motivated models (DQ and 
"heavy-bag") of the quark phase matched with real- 
istic RMF models (either "hadron-MWl" or "hadron- 
MW2") of the hadronic phase predict the onset of the 
phase transition at very high critical densities (typi- 
cally n h c ~ 10n ). Note that both "hadron-MWl" and 
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FIG. 4: (Color online) Effective quark masses (left panel) and bag parameters (right panel) for the DQ models (curves), cf. 
Eq. JKJ, and bag models (horizontal lines) as functions of n/no for T = at /3 equilibrium. Notation is the same as in Fig. 




FIG. 5: (Color online) Pressure as a function of the chemical potential (left panel) and n/no (right panel) for T = at j3 
equilibrium. Notation for the quark models is the same as in Fig.0 Band boundaries for the "heavy-bag" and "light-bag-155" 
model correspond to c = and c = 0.3: from top to bottom in left panel and from bottom to top in right panel. Hadronic EoS 
of the conventional Walecka model ( "hadron W" ) is displayed by long-dashed line. Two versions of the hadronic model of Ref. 
[5^ |. i.e. "hadron-MWl" and "hadron-MW2" , are presented by dashed-dotted-dotted and dotted lines, respectively. 



"hadron-MW2" models do not permit neutron stars with 
central baryon densities > 8 no, see Ref. Thus, 
our lattice- Q CD-motivated models (DQ and "heavy-bag") 
matched with realistic RMF models ( either "hadron- 
MWl" or "hadron- MW2") do not allow the existence 
of hybrid stars. The stiffer the hadron EoS, the lower 
the critical density is. However, the central den- 
sity of the limiting mass star also decreases. In the 
case of the DQ and "heavy-bag" models matched with 
the "hadron- W" model representing unrealistically stiff 
EoS (i.e the stiffest EoS among those still used in cal- 
culations), the typical critical density w Auq. The 
critical masses M(n cent = n'*) for the occurrence of 
the hybrid stars are listed in Table IIVI Three hadron- 
quark model combinations, i.e. "hadron-W" + "l-loop" 
and "hadron- W" + "heavy-bag" with c = and c = 0.3, 
predict the existence of hybrid stars with M > 2.7 M©. 



These mass values are only slightly below the limiting 
mass. The hadron-quark combination "hadron- W" + "2- 
loop" predicts no hybrid stars. 



The "light-bag-155(c = 0.3)" and "light-bag-200(c = 
0)" models result in critical densities in a range which 
is usually considered in papers devoted to hybrid stars. 
For the "light-bag-155(c = 0)" model matched with any 
of the considered hadron models, quark matter is always 
more stable than hadronic matter (the quark pressure is 
higher), and hence quark stars of an arbitrary size are 
possible. In this case only a softer (compared to those 
we have used) hadronic EoS could allow the possibility 
of hybrid star existence. 
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hadr.\ quark 


"1-loop" 


"2-loop" 


heavy-bag(c = 0) 


heavy-bag(0.3) 


light-bag-155(0.3) 


light-bag-200(0) 


W 


4.4/20.3 


4.8/19.3 


3.9/12.1 


4.4/12.6 


1.3/2.6 


2.2/6.4 


MW1 


9.8/27.7 


11.2/29.7 


9.3/22.0 


10.8/24.1 


3.7/4.3 


4.4/7.1 


MW2 


12.2/31.3 


14.1/34.5 


11.5/25.1 


13.7/29.1 


5.8/7.1 


5.4/9.1 


TABLE III: 


The hadronic and quark n\ critical densities (in units no) according to the Maxwell construction. 


hadr.\ quark 


"1-loop" 


"2-loop" 


heavy-bag(c = 0) 


heavy-bag(0.3) 


light-bag-155(0.3) 


light-bag-200(0) 


W 


2.7 




2.8 


2.7 


0.6 


2.1 


MW1 










1.6 


1.8 


MW2 










1.8 


1.7 



TABLE IV: Masses of neutron stars M (in units Mq) corresponding to the critical densities in the hadron phase, , listed in 
Table EH 



C. Quark stars 

Now we turn to the question of the possible existence 
of pure quark stars. As argued in Ref. the quark 
star may be surrounded by only a tiny crust separated 
from the core by a strong electric field. Therefore, we will 
neglect contribution of the crust. Then the star mass and 
its radius result from the Tolman-Oppcnheimer-Volkoff 
equation (in the standard notation, cf. 



dP 

dr 
dm 
dr 



Gem 1 + P/e 



r 2 1-2G m/r 



AnPr- 



(25) 



4wr e 



which is solved for a given e(P) dependence with ini- 
tial conditions m(0) = and P(0) = P C cnt, where P C ent 
is a central pressure corresponding to a central density 
?icentr- The latter is the input parameter of the calcu- 
lation. For the sake of simplicity, we neglect effects of 
possible rotation of the star, which however may affect 
the mass-radius relation of the newly born stars and also 
influence the dynamics of the possible phase transition of 
the neutron star to the hybrid or the quark star, cf. [H^ . 

The pressure as a function of the energy density is 
shown in Fig. for the models under consideration. We 
see that for all the quark models these dependencies 
can be well parameterized with appropriate accuracy by 
straight lines e = aP + b, with b = AB, where B has the 
meaning of an effective bag parameter, cf. Ref. |2lL |22|. 
For the "massless-bag" model, i.e. with all quarks being 
massless, we would obtain a — 3(1 — c), b = 4B. Values 
of the effective bag parameter B and slop-coefficients "a" 
are listed in Table [V] Note that the quasiparticle model 
of Refs HI E2 resulted in 3.1 < a < 4.5 and B 1 ' 4 > 
200 MeV (typically B 1 ' 4 ~ (230^-240) MeV). Within our 
lattice-QCD-motivated models (DQ and "heavy-bag") 
we obtained essentially higher values: a = 5 -j- 6 for and 
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FIG. 6: (Color online) Pressure versus energy density at T = 
in /3-equilibrium for various EoS's. Band boundaries for 
the "heavy-bag" model correspond c = and c = 0.3 (from 
bottom to top). Notation is the same as in Figs Q and |5] 



B 1 / 4 = 212^-304 MeV. As seen, the e(P) lines are quite 
different for the "light-bag" , "heavy-bag" and DQ mod- 
els. Due to that one may expect quite different gravita- 
tional properties of objects described by these models. 

In Fig. [7] we show mass-radius (M-R) relation for our 
models of quark stars, M = m(r = R). Two extreme 
cases are demonstrated: quark matter in the normal 
phase (left panel) and the CFL phase of quark matter 
with a high pairing gap A = 100 MeV (right panel). In 
the latter case the contribution of pairing to the energy 
density is the largest. For other superconducting phases, 
the results will range in between these two extremes pre- 
sented here. Note that occurence of high n q c values in 
our DQ and "heavy-bag" models, see Table ITTT1 implies 
that the CFL phase takes place in these cases. Fig. |S] 
demonstrates the quark star mass as a function of the 
central baryon density, again for the normal phase with 
A = (left panel) and the CFL phase with A = 100 MeV 
(right panel). The limiting masses and central densities 
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Model "1-loop" 


"2-loop" 


heavy-bag(c = 0) 


heavy-bag(0.3) 


light-bag-155(0) 


light-bag-155(0.3) 


light-bag-200(0) 


a 5.01 


5.97 


5.89 


5.41 


3.11 


3.09 


3.15 


B 1/4 , MeV 304 


293 


216 


212 


159 


159 


202 



TABLE V: Coefficients of the linear interpolation of the e(P) dependence. 
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FIG. 7: (Color online) Quark star mass - radius relation of pure quark stars for the normal quark phase with A = (left panel) 
and the CFL quark phase with A = 100 MeV (right panel). Notation is the same as in previous figures. Band boundaries c = 
and c = 0.3, from bottom to top. 



are listed in Table IVII for normal quark matter and in 
Table EH for the CFL quark matter. Note that CFL 
strangelets were discussed in Ref. [^J • 

In the case of "light-bag-155" model the pairing may 
substantially affect the EoS and hence quark star masses. 
The difference between normal and CFL quark star 
masses increases up to 0.6 M Q in the region M ~ Mu m . 
The limiting mass is Mu m ~ (1.6-5-2. 3)Af© and the radius 
is ~ 9-5-12 km. Thus, for these models the M-R and M- 
Hcent relations are similar to those one usually uses for 
pure hadronic EoS's, in agreement with the statement 
of Ref. For the "light-bag-200(c = 0)" model the 

limiting mass is Mu m ~ (0.9 -5- 1.2) M© and the radius 
is ~ 6 km. This limiting mass is smaller than several 
well measured masses of pulsars. Thus we may conclude 
that the "light-bag-200" model does not describe ordi- 
nary pulsars but rather another family of compact ob- 
jects. The effect of the pairing is less pronounced here 
than for the "light-bag-155" model. Let us remind that 
the "light-bag-155(c = 0.3)" and "light-bag-200 (c = 0)" 
models matched with either "hadron MW1" or "hadron 
MW2" hadronic models also predict existence of hybrid 
stars. Had we assumed a softer hadron EoS, the corre- 
sponding hybrid stars would become impossible and we 
would only have pure quark stars. 

For the "heavy-bag" model and especially for the DQ 
model, the masses and radii are still smaller and the in- 
terior region, denser. For the "heavy-bag" model the 
limiting mass is M ~ (0.6 -5- 1) M© and the radius is 



(~ 4 -5- 6) km. These configurations are similar to those 
previously studied in the model of Ref. 22]. For the 
DQ model the limiting mass is M ~ (0.4 ~ 0.5) M© 
and the radius is (~ 2.5 -5- 3) km. The central density 
is ~ (50 -5- 60)no. Such dense quark matter is probably 
realized in the CFL phase that has peculiar cooling prop- 
erties. Small and dense objects might be of interest with 
respect to MACHO events [H^. 



V. CONCLUSIONS 

Models for deconfined quark matter are extensively 
used in astrophysical applications. The discussion of hy- 
pothetical hybrid and quark stars is a rather hot topic. 
In this paper we emphasize that it is essential for a de- 
confined quark model to be constrained by the lattice 
QCD results, which became recently available at finite 
baryon chemical potential. Rapidly progressing lattice 
calculations may essentially restrict models used for the 
description of cold and dense quark matter. In this pa- 
per we are trying to make a step in this direction. At 
this time, we would like to draw attention to problems 
with traditionally used simple standard bag and, prob- 
ably, NJL models. They are not able to reproduce the 
lattice data without significant modifications. 

We have considered the dynamical quasiparticle 
model [2^] and the heavy-bag model, which fit the lat- 
tice results. The heavy-bag model uses heavy masses 



12 




FIG. 8: (Color online) Quark star mass - central density relation. Notation is the same as in Fig. Q 



Model 


"l-loop" 


"2-loop" 


heavy-bag(c = 0) 


heavy-bag(0.3) light-bag-155(0) 


light-bag-155(0.3) 


light-bag-200(0) 


M Hm 


0.43 


0.43 


0.62 


0.67 1.60 


1.62 


0.99 


Tkcentr 


61.0 


55.6 


28.3 


24.3 9.42 


8.54 


19.8 



TABLE VI: Limiting quark star masses (in units Mq) and corresponding central densities (in units no) for the normal phase 
(A = 0) of quark matter. 



of quarks and gluons to reproduce the lattice results 
well. We have demonstrated, that these models, being 
matched with realistic hadron models of EoS, rule out 
the existence of hybrid stars. This implies that families 
of ordinary neutron stars and pure quark stars are well 
separated. Within these models the quark stars have low 
masses (M < 1 Mq), are dense (n > 10no) and very 
compact (R < 6 km) . This result complies with that ob- 
tained in Ref. by means of a different extrapolation 
of lattice data to cold dense baryonic matter. If such 
a quark star existed in a binary system (together with 
a white dwarf, a neutron star, or together with another 
quark star), the orbiting would be different from that ex- 
perimentally known for both white dwarf-neutron star 
and double neutron star systems. Moreover, such dense 
quark matter may be in the color-flavor-locked supercon- 
ducting phase and hence possess peculiar cooling prop- 
erties. Other properties of these pure quark stars also 
differ from those of neutron stars that may give a chance 
to distinguish them experimentally, e.g., see Ref. [57I . 

Another interesting question arises, what could be a 
possible formation mechanism of such compact and dense 
quark stars? For the birth of a quark star, a hybrid 
star configuration is first required in order to produce a 
quark core either during supernova collapse or during the 
long-term accretion of the matter in a binary system, or 
in the course of the merging of two neutron stars. The 
nucleation dynamics of quark matter in a cold neutron 
star was considered recently in Ref. (5^- Then also a 
mechanism to blow off the hadron shell is needed, e.g., 



like that discussed in Ref. 60] for a strong first-order 
phase transition. 

In the present paper the models for deconfined mat- 
ter, which are fitted to lattice predictions, allow neither 
hybrid stars nor quark self-bound objects (strangelets). 
We have demonstrated that only the unrealistically stiff 
original Walecka model ( "hadron- W") matched with ei- 
ther the DQ or "heavy-bag" models may allow for a 
small quark core. Thus, we actually see no appropri- 
ate mechanism to form quark stars that are described by 
our DQ and heavy-bag models. Temperature, rotation, 
and neutrino trapping effects should be included. They 
may stiffen the hadron EoS |6l| and perhaps may help 
to resolve the problem. Some exotic ways to produce 
quark stars through either evolution of white dwarfs or 
inflation stage at the quark-hadron phase transition are 
discussed in Refs [53, |62j . Certainly, this problem still 
needs further consideration. 

It is necessarily to mention uncertainties in the present 
considerations and the conclusions derived. First of all, 
they are connected with lattice data 14], which we ex- 
trapolate to the domain of cold and dense baryonic mat- 
ter. These data obtained on a discrete lattice are only 
poorly extrapolated to the continuum limit. Another 
problem with the lattice data consists in the poor repro- 
duction of the chiral limit. As we have mentioned already, 
the pion turns out to be very heavy in lattice calculations, 
e.g. ~ 770 MeV (3^. All this may result in certain 
errors in the parameters of our effective models, deduced 
from these data. We have tried to estimate this kind of 
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Model 


"l-loop" 


"2-loop" 


heavy-bag(c = 0) 


heavy-bag(0.3) 


light-bag-155(0) 


light-bag-155(0.3) 


light-bag-200(0) 




0.47 


0.49 


0.85 


0.96 


2.13 


2.26 


1.17 


T^centr 


53.2 


46.9 


16.1 


13.0 


6.41 


5.02 


15.8 



TABLE VII: The same as Table EH but for the CFL phase of quark matter with A = 100 MeV. 



uncertainties by applying different extrapolations from 
the lattice data, i.e. the "l-loop" and "2-loop" versions 
of the DQ model and the "heavy-bag" model. Indeed, 
whereas the "2-loop" version closely reproduces the lat- 
tice data the "l-loop" version slightly (by 10%) mis- 
fits its normalization, thus leaving room for uncertainty 
in the extrapolation to the continuum and chiral limits. 
As for the "heavy-bag" model, it covers the lattice data 
by a rather broad band (due to variation of c) . 

On the other hand, the grounds for the standard bag 
model, in our case represented by the "light-bag-155" and 
"light-bag-200" models, are not entirely reliable. Indeed, 
the bag parameters of standard bag models are fitted to 
reproduce the masses of hadrons, i.e. few-body systems, 
and then applied to describe matter. It could happen 
that the situation here is similar to that of Walecka model 
for nuclear matter. If we fit the Walecka model to repro- 
duce the deutron or helium and then apply it to 208 Pb, 
we surely get nonsense. Therefore, predictions of these 
models for the EoS should be treated with great care. 

Theoretical uncertainties in our predictions should be 
mentioned as well. The region of the vicinity of Tq de- 
serves special discussion. Besides the above mentioned 
uncertainties of the lattice calculations, as it is known, in 
the (2+l)-fiavor case the phase transition at small chem- 
ical potentials is of the cross-over type . This implies 
that the hadronic and quark-gluon degrees of freedom es- 
sentially interact in this region, e.g., see Ref. [63^, This 
is why we cannot pretend to reproduce the lattice data 
in this region by means of only quark-gluon degrees of 
freedom incorporated in the DQ model. Our aim was 
to fit these data beyond this vicinity region and then to 
extrapolate them to the domain of the cold and dense 
baryon matter. As it is commonly expected, the phase 
transition in this domain is of the first order, which im- 
plies a less severe role of the hadron-quark-gluon inter- 
action. Therefore, we hope that the method of matching 



hadronic and quark-gluon EoS's, applied in this domain, 
is reasonable. However, the uncertainty of such an ex- 
trapolation should be kept in mind. We cannot also ex- 
clude that our models are oversimplified and some impor- 
tant physics is left aside, e.g., the model suffers of a lack 
of the chiral symmetry restoration, thermal fluctuations 
are not included, which may play an important role in the 
vicinity of the critical temperature but die out at smaller 
and higher temperatures, etc. This may somehow modify 
the fitting parameters of models. Quantum fluctuations 
like instantons and gluon condensate may partially sur- 
vive and additionally contribute in the low temperature 
region, cf. Ref. j64j. Thus, it might well be that the 
DQ model, as well as the heavy-bag model, should be 
replaced in the future by more sophisticated and realis- 
tic microscopically based models. Nevertheless, a certain 
confidence in our results is supported by the fact that 
other independent extrapolations of the lattice data to 
the cold dense baryonic matter [22Ll24| resulted in rather 
similar conclusions. However, it would be of interest to 
search for new theoretical approaches for extrapolating 
the EoS from T > 190 MeV to T = and from = 530 
MeV upwards. 
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